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spreading control sea-water intrusion, the author states: (paragraph following 
Eqs. 17) maintain control indefinitely the rate infiltration must such 
cause rise the mound least mean sea For unconfined 
homogeneous isotropic aquifer the mound elevation must exceed that mean 
sea level develop effective equilibrium condition repelling sea-water 
intrusion. According the Ghyben-Herzberg the greater speci- 
fic gravity salt water will balance correspondingly greater head fresh 
water saturated ground-water zones. Following the author’s notation, this 
may expressed 


which equals the specific gravity fresh water; the gravity 
salt water; and AH, the incremental head fresh water required the mound 
maintain equilibrium. Solving for AH, 


Qs — OF 


illustrate the magnitude the increased fresh-water head required, 
inserted Eq. 126, together with 1.026 and 1.000. 
value 3.5 Thus, assuming that undiluted sea water 
direct contact with the fresh ground water, the fresh-water mound must 
maintained more than above mean sea level control sea-water intrusion. 

The author has limited his analytical treatment and model studies cases 
unconfined ground-water conditions only, and has not investigated the problem 
pressure aquifers. However, the impression given that important 
application spreading that control sea-water intrusion. This applica- 
tion important coastal areas where unconfined aquifers exist; but 
California where the problem sea-water intrusion becoming increasingly 
serious, most the contaminated basins are located pressure areas. 
Table are the known localities sea-water intrusion California 
1950; the areas acres underlain salt water; and the type aquifer—con- 

paper Paul Baumahn was published August, 1951, Proceedings-Separate No. 86. 


The numbering footnotes, illustrations, tables, and equations this Separate continuation the 
consecutive numbering used the original paper. 
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fined unconfined. Seven the known areas sea-water intrusion are 
situated pressure aquifers, three non-pressure aquifers, one semi-con- 
fined aquifer, and two undetermined aquifers. the 90,525 acres that have 
been found underlain sea water, 87,965 acres occur pressure areas. 
control intrusion pressure areas, methods other than spreading must 
employed. Current studies include the technical and economic 
feasibility constructing artificial subsurface barriers, and also creating 
feash-water pressure ridges above mean sea level means injection wells. 
The evaluation the mound width developed circular spreading ground 
above inclined impervious stratum (Eq. 124) showed that circle 500-ft 
radius provided mound width approaching 10,000 ft. The same equation 


TABLE AREAS AND AQUIFERS INVADED 


Ground-water basin valley (acres) Type aquifer Authority 

Central Coastal Plain............ 3,500 Confined Zander and 
West Coastal 16,000 Confined Zander and 
Sacramento-San 2,400 Semi-confined Tolman* 

Total known area............. 90,525 


Area underlain salt acres. Water Supply Paper No. 446, U.S. Geological Survey, 1919. 
Bulletin No. Div. Water Resources, State California, Approximately 500 Bulle. 
tin No. 16, Div. Water Resources, 1933, Bulletin No. 52, Water Resources, State California, 
1946. Between Salinas Valley and Pajaro Valley. “Santa Cruz-Monterez Counties Investigation,” 
Bulletin No. State Water Resources Board, State California (publication pending). Water Supply 
Paper No. 519, Geological Survey, 1924. Joaquin Valley between Pittsburgh and 
Tolman, McGraw-Hill Book Co., Inc., New York, Y., 1937, pp. 


378-379. 


may applied determine the mound width resulting from recharge well. 
This would have application consideration were being given the main- 
tenance fresh-water barrier along coast unconfined aquifer using 
line recharge wells parallel the coast. mound width would determine 
the spacing the wells. Since this width independent the spreading area 
size and depends only the rate aquifer recharge, and the physical 
parameters and estimate can made based representative rate 
well recharge. Two published accounts giving quantitative measurements 
recharge wells southern show average recharge rates varying 
from 0.52 2.92 acre-ft Taking value acre-ft per day 
example and using the author’s values 100 per hr, and 


Investigational Work for Control and Prevention Sea-Water Intrusion into Ground 
Water Banks and Bookman, Report State Water Resources Board, Div. Water 
Resources, State California, 1951. 

Water for Storage Underground,” Mitchelson and Muckel, Technical 
Bulletin 578, Dept. Agri., Washington, C., 193 

Lane, Transactions, Am. Geophysical Union, Vol. 15, 1934, pp. 523-527. 
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0.01, Eq. 124 yields mound width, 54.5ft. Actual wells would need 
spaced even closer than this distance develop suitable mound height 
and leave intervals between wells without fresh-water flow. this 
example, the mound width would have reached only 54.0 distance 
1740 downstream from each well. The spacing suggests that economic 
considerations alone would govern applying this method control sea-water 
intrusion. Having line recharge wells with only some between each 
well indicates that recharge spreading would better method actual 
situation approximating these conditions. 


Max ASCE.—This paper excellent mathematical treat- 
ment the problem water recharge, problem that becoming more and 
more important and which the present literature relatively scarce. The 
author shows how exceedingly complicated the mathematical part the purely 
hydraulic phases can become skilled hands unafraid hard work. 
practical applications ground-water recharge many additional factors appear, 
such nonuniformity the soil and basin boundaries, variations the 
chemical and physical character the recharge water, silting the soil (either 
silt the recharge water effect the soil particles), pollution 
ground water, influences the temperature and chemical composition the 
ground water, predictions rate inflow, duration inflow and periods 
between recharges, pumpage the vicinity, and backwater effects and general 
variations ground-water levels due rainfall and floods. Many these 
factors cannot represented mathematical formulas they depend too 
much local conditions. good, therefore, have some calculating meth- 
ods for those phases that are amenable mathematical treatment and this 
respect this paper marks great step forward. 

Much the mathematics used the paper probably unfamiliar many 
practical engineers, few them had remember differential equations. 
may asked, therefore, all this higher mathematics necessary for the 
solution the problem, and justified view the simplified assumptions 
that are made. the following, only few examples will discussed show 
the general trend answer the foregoing question. treat every devel- 
opment the paper the same manner would fill book. 

Eq. 10a checks the writer’s Fourier series expansion provided the member 
dropped the general formula for the Fourier series— 


n=1 


The use the Fourier series expansion makes Eq. look formidable when 
actually nothing else than 


obtain this result not necessary through the Fourier series ex- 
pansion. 
Head, Eng. Research, State Water Survey, Urbana, 


th 
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Applying condition Eq. 


and, this special value introduced into Eq. the result derivation 
Eq. 128 without the detour the Fourier series expansion. The develop- 
ment Eq. 128 based Eq. the approximate equation the stable 
mound att Referring Fig. Eq. the actual formula for y’, not 
somewhat simpler expression than Eq. 128 obtained using Eq. 
replaced 


Similarly, the place Eq. 14, 


the same procedure used the derivation Eq. 128 applied Eq. 130, 
replacement Eq. 128 becomes 


The principles here developed can expanded other formulas the 
paper, and can seen that, without recourse the Fourier series expansions, 
formulas are obtained that are simpler form and easier for numerical 
tions than those the paper. This important the Fourier series, although 
beautiful theory, are notoriously cumbersome for numerical calculations. 

may also questioned the straight-line slope the correct approxima- 
tion for the steady mound. based the assumption Jules Dupuit 
that the velocity vertical 
and amount equal 
that defined the Darcy 
for flat slopes for the 
sake integrality the 
equations. 

The Dupuit assumptions 
are nonpotential but give 
close approximation flat 
slopes, where the sine the 
slope can replaced the 
tangent. The following de- 
velopments are based these 
assumptions and show that even these simplifications furnish curves the gen- 
eral character obtained the potential theory and Mr. Baumann’s experi- 
ments. Assuming point source with circular runoff around it, the areas are 
follows: 


— 
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distance 


Therefore, the change area 
—neglecting the higher differential (see Fig. 22). 


Inasmuch 


TABLE THREE CASES WATER SURFACE CURVES 


Case Value Area Velocity Slope Water surface 
Constant Constant Constant Straight line 
Positive Increasing Decreasing Flattening Concave line 
Negative Decreasing Steepening Concave line 


there are three cases for the water surface curve (see Table 2). When these 
cases occur? Since 


Integrated, 
and 
(139) 


For flat slopes but the actual value not yet known. Under the 
conditions assumed the paper 


(140) 
Since— 
and 
(142b) 


a 
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The condition exists therefore only the immediate vicinity 
elevation fixed natural conditions. The value that point can 
determined study the curve the water surface. Integrating Eq. 
141, equation the steady mound derived: 


For infinitely large spreading ground, and for limited spreading 


ground log these two cases, Eq. yields 


and 


The point with therefore also point inflection and its 
coordinate, from the condition defined Eq. 133, computed 


Therefore, this point lies the line and 138 


This steep slope explains why, actual field and laboratory observations, 
there sharp drop near the recharge area. 

slope 45° never observed because, the field and the model, there 
are influences that flatten the slope—namely 


(a) For steep slopes Dupuit’s assumptions not yield good approxima- 
tion; and 

(b) There area inflow instead point recharge. Eq. 146 shows 
two additional properties the surface configuration 
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such curve convex between the source recharge and point the 


other hand, 


and such curve concave from point outward. Considering these proper- 
ties the curve the steady mound does not seem appropriate apply 
refined computation involving assumption that the water surface can rep- 
resented straight line. 

The author presents curves the type here developed, but also inverse 
curves. However, the writer’s belief that the drop shown near the outflow 
due backwater less height than the critical 


number cases when successful structures were built and engineering oper- 
ations performed from experience and partly intuition, but without adequate 
theoretical background the part engineers. Theoretical clarifications 
based physics and were given generally much later scientists 
and scientifically inclined engineers other than the original ones. Such, for 
instance, was the case ancient earth dams that were built since times im- 
memorial without their builders having any notion soil sub- 
surface flow net theories. Conversely, the case the present paper have 
extraordinary example modern engineer who for many years was engaged 
the technical direction spreading operations large scale Southern 
California and then himself produced complete theory spreading tending 
explain the minute details the phenomenon. Under given circumstances 
exceedingly difficult discuss this interesting paper, not only because 
the abundant and rather complicated mathematical analysis presented, but 
chiefy because the fact that the author, virtue his practice, strongly 
prepared defend his views. The latter, addition, are backed well-per- 
formed experiments which the writer had opportunity examine personally 
1948. The writer, therefore, wishes limit his discussion mostly the 
comparison the author’s philosophy approach the given problem with 
his own views. 

After having read Mr. Baumann’s paper the writer became convinced 
that least three problems the underground flow belong the same class. 
These are—the theory consolidation; drainage the pavement base 
and spreading. From these three the spreading problem the most compli- 
cated because this case the flow confused especially the experimental 
stage, capillarity, mostly above the saturation line, and also the boundary 
conditions, the bottom the experimental box. fact, the latter condi- 
tions are indicated the travel fringes marked and on, minutes 


Charles Jaeger, Basel, Switzerland, 1949, 366. 


Cons. Engr., San Francisco, Calif. 
Course Drainage for Airport Casagrande and Shannon, Proceedings- 
Separate No. 75, ASCE, June, 1951. 
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Figs. and 20. the opinion the writer these fringes are produced 
the colored liquid that falls the bottom the box through water and spreads 
horizontally due the combined action hydraulic head and momentum. 
The writer observed somewhat similar phenomenon his experiments 
horizontal cylindrical tubes filled with dry sand. the latter 
case water just under small head first very quickly spreads horizontally 
along the bottom the tube and afterwards invades the tube capillary 
way moving The writer believes that the use parabolic bottom 
the experimental box (Fig. 21) was fruitful idea since such way the 
true character flow nature hetter duplicated the experiment than 
the case horizontal bottom the experimental box. 

The solution the three aforementioned problems based one way 
another the well-known differential equation heat transfer through 
prismatic, nonradiating bar with constant coefficient proportional the 
coefficient permeability the case subsurface flow the coefficient 
not constant, however, and gradually decreases either the process 
drainage spreading advances. This pronounced phenomenon 
the case soils with fine-grained admixtures. The situation better the 
case flow through coarse pervious materials since this case relative 
constancy the coefficient permeability during drainage spreading may 
expected. Hence, the author’s statement that spreading grounds 
are located soils that approach idealized conditions extending over great 
encouraging. noticed that the same favorable condi- 
tion prevails the case runway bases made gravel similar materials. 

The hydraulics the open-channel flow had existed before the advent 
the theories the subsurface flow; and undoubtedly the latter have been 
influenced the former. Thus, the theories quasi-horizontal subsurface 
flow the existence horizontal slightly inclined impervious bottom (like 
those occurring some rivers) predicated. reality, however, the earth 
mass may pervious indefinite depth and saturated with ground water. 
The most realistic explanation why impervious bottom introduced the 
theories quasi-horizontal subsurface flow the necessity have plane 
from which the ordinates the subsurface stream could measured; other- 
wise the formulas cannot derived. Such the case the Dupuit formula, 
the theory the pumping test, the theory the drainage pavement 
and, finally, the theory spreading presented this paper. 

The writer feels, however, that, the case very deep ground-water 
stream with large values spreading would affect the lower part this 
stream very little, all. stated sometimes soil that that 
water has shearing strength, but this statement relative value only. 
The very concept viscosity based the possibility that shearing stresses 
exist viscous fluid such water. The hydraulic gradient pushing the 
spreading stream down larger than that the underlying ground water; 
hence, the spreading stream moves with larger velocities than the underlying 
ground water. This causes shearing stresses and corresponding increase 


Capillary Phenomena Sandy Krynine, Proceedings, 29th Annual 
Meeting, Highway Research Board, Washington, C., 1950, pp. 520-530. 
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velocities the ground water, which fact permits gradual incorporation the 
mound wave into the original ground-water stream. The shearing stresses 
have decrease rapidly from the surface the original ground water down- 
ward and depth become negligible for all practical 
pattern may easily visualized analogy from the diagram shearing 
stresses body with horizontal slightly inclined upper surface over 
which heavy object with plane bottom surface sliding. The writer’s 
reasoning valid, his opinion, least for the first and second phases the 
spreading process. 

conclusion the writer wishes congratulate Mr. Baumann his as- 
siduous and conscientious work which represents valuable contribution the 
mathematical theory spreading. 


several years, the writer had principal objectives, namely—to find sound, 
theoretical interpretation observed ground-water phenomena related 
spreading, and clarify attendant, fundamental concepts which appeared 
afflicted, some degree least, confusion. Above all, was felt that 
the wave phenomenon and the significance controls had not been fully 
evaluated and recognized, respectively. 

Mr. Todd seems interpret the writer’s expression, least mean sea 
falling short possibly 3.5 above mean sea level. Such was not 
the intention. Since 1950 the writer has been charge investigation, 
the Los Angeles County Flood Control District, means for the 
stemming sea-water intrusion the West Coastal Basin. This basin (which 
Mr. Todd refers Table Coastal has shore line about 
miles length along which about two thirds the principal aquifer, known 
Silverado Zone, are confined and the balance unconfined. The latter 
predominantly located the southern part the area. 1950 spreading 
tests were conducted confirm the unconfined character the aquifer there, 
determine the permeability the soil, and ascertain the salinity the 
ground water. The investigation led the conclusion that, this area, sea- 
water intrusion could best controlled means liquid barrier, created 
through spreading from strip area from shallow gravel-filled injection pits, 
apart. This barrier would have located well inland from the sea- 
shore have the assurance control. Therein lies distinct difference 
between the mechanics liquid barrier open and confined aquifer. 
The latter, Mr. Todd correctly states, not treated the paper, although 
recharge test was conducted that same year abandoned well the con- 
fined zone the aquifer. Both tests plus third one the northern quasi- 
unconfined area were described report the Los Angeles County Flood 
Control District, prepared the spring 1951, together with study the 
location and components complete liquid barrier over the entire 12-mile 
length. Finally, included recommendation for pilot test recharge 
wells, extending over front the confined area, some 2,000 inland 
from the seashore. Although the chloride concentration this location may 


Asst. Chf. Engr., Los Angeles County Flood Control Dist., Los Angeles, Calif. 


BAUMANN GROUND-WATER SPREADING 


correspond coefficient 1.01, and may call for rise the barrier above 
mean sea level the chloride concentration the proposed injection 
pits the open aquifer very much smaller, require rise there only 
few tenths foot above mean sea level. the wording the paper 
was meant imply that, whereas some locations the rise the mound above 
mean sea level would negligible stem intrusion, could appreciable 
other locations. Where the latter happens the case, each foot extra 
rise may indeed represent considerable part the effort. Hence, Mr. 
Todd’s point well taken. 

Fortunately, the conditions the open area the West Coastal Basin are 
not difficult they would the physical parameters used the writer’s 
prototype example applied. These values actually apply spreading ground 
debris cone area and not the open aquifer the West Coastal Basin. 
The respective values would more nearly: 0.25 per 900 
ing mound width, ft, would assure complete closure based the 
aforementioned pit spacing ft. 

reading Mr. Suter’s able discussion, one should bear mind two funda- 
mental premises: (1) The writer confined linearization the stable station- 
ary mound the two-dimensional case only and (2) “the improbability the 
existence circular control” was emphasized connection with the three- 
dimensional case. The linearization the two-dimensional mound per- 
missible approximation long accordance with the basic premise. 
The discrepancy between the linearized and the true, stable mound shown 
Fig. exaggerated because the comprehensive model tests the afore- 
mentioned, basic premise was not satisfied. fact, the height the stable 
mound was nearly one half the initial ground-water depth. This exaggera- 
tion was intentional facilitate observation the one hand and ascertain 
the effect this deviation from the basic premise the theoretical results. 
well-known fact that mathematics often funishes valid results well beyond 
the imposed limitations. 

Mr. Suter makes commendable effort simplification Eq. 
introducing boundary condition Eq. Apparently the statement, 
still does not satisfy condition was thereby overlooked. 
Eq. does not satisfy condition then condition does not satisfy Eq. 
Hence, necessary, first all, transform Eq. into equation which does 
satisfy This was accomplished through regrettable, 
course, that condition singularly difficult satisfy and that the elegant 
operation suggested Mr. Suter does not answer the purpose. Unfortunately, 
leads the result expressed through Eq. 128, which describes family 
straight lines, radiating from the control toward the spreading ground. Hence, 
accordance therewith, not only are the undisturbed ground-water surface 
and the surface the stable mound linear but also all the intermediate surface 
configurations the growing mound. Again, referring Fig. 14, obvious 
that such not the case. All equations derived from Eq. must necessarily 
satisfy the latter. This true all the writer’s equations for the growing 
and the disappearing mound, but not true Eq. 128. 
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and which express the entire, two-dimensional phenomenon are 


not formidable they may appear first glance. La, and 
Eqs. check Eq. Length, La, may any finite distance. Hence, the 
determine the progress and shape the wave both the upstream and down- 
stream direction for such arbitrary time elements correspond the first 
phase. This shown Fig. which based the same numerical values 
the dash lines Fig. except for instead ft. shows that 
this case the time required for the wave reach the control, and for discharge 


Ld w > 


Height Mound, Inches 


Initial Ground 
Water Depth 
Inches 
(Not Scale) 


Distance, Feet 


addition the initial ground-water flow commence, the order 
min compared 0.5 min the first case. also shows the flexibility 
equations this type. For Eqs. also check the analogue 
the for heat flow nonradiating, prismatic bar semi-infinite 
length which reads (in terms ground-water parameters) 


—z2 7-1/4 
= — DP, 


v o 


Ingersoll, Zobel, and Ingersoll, McGraw-Hill Book Co., Inc., New York, 
Y., 1948, 110, 
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integral which readily determined from Eq. 150 reduces 


which also checks the corresponding term derived from Eqs. 17. 

Mr. Suter states correctly that the writer’s equations for are actually those 
for according Fig. and that the relations shown 130 and 131 
should used for the respective transformation. Unfortunately, however, this 
transformation does not free Eq. 132 the discrepancy attending Eq. 128. 

the light the foregoing, Mr. Suter’s check the validity the 
writer’s linearization does not seem pertinent that his development Eqs. 


133 and including 148, entirely based three-dimensional flow for 


Actually the writer’s formulas Part II, Section correspond the afore- 
mentioned ones, except that Mr. Suter assumes stationary extracted, flow 
(positive right side Eq. 141) and therefore circular control, whereas the 
writer’s development based unsteady injected flow (negative right side, 
Eq. 107) and control. Furthermore, Mr. Suter assumes point source 
(Fig. 22) instead circular spreading ground used the writer 
(Fig. 9). once apparent that reduction the radius the spreading 
ground such approach point source would result surface configura- 
tion quite different from that shown Fig. 22. 

setting the equations for the two cylindrical areas distance 
apart, preceding Eq. 133, Mr. Suter introduces negative increment for 
positive increment dz, whereas both should entered positive. 
mathematics tell the end whether positive negative. 

Hence, Eq. 138 for the condition should read 


which defines family square hyperbolas with the y-axis and line parallel 
to, and distance below, the asymptotes. The intersection any 
these hyperbolas with the undisturbed ground-water surface corresponds 
the distance which the circular control assumed. This solution satisfies 
the fundamental law that impossible maintain all three—area, velocity, 
and slope constant—for increasing These hyperbolas are special case 
the so-called ‘‘cone important emphasize, however, 
that recharge well cannot treated water well is, 
assuming stationary flow Mr. Suter did. The reason that water well 
itself control whereas the recharge well not. Hence, the latter case 
the flow unsteady and radiates from the well uniformly any distance 
form wave without ever becoming stationary within finite limitations. 
the course this performance the rate recharge must diminish illustrated 
Fig. (dash line). This diminution the rate recharge is, purely and 
simply, necessity; yet, has been attributed causes 


Ingersoll, Zobel, and Ingersoll, McGraw-Hill Book Co., Inc., New York, 
Y., 1948, see appendix. 
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other than the real one, recent investigators. The writer means 
oblivious the fact that causes entirely different nature exist; but 
believes taking the first step first and taking the probable into considera- 
tion ahead the possible. 

Mr. Suter’s final paragraph, reference made the critical depth the 
control accordance with Charles Jaeger’s which compares 
with Eq. equations are approximations. Contrary the writer’s 
development Section Part Mr. Jaeger arrives his value through sub- 
stitution circular ares for the parabolic equipotential lines. This leads 
sharp drawdown the phreatic line and vertical tangent the control. The 
point contact this tangent determines the critical depth. either case 
the critical depth defined the minimum height the phreatic line the 
control regardless how far below the backwater surface may be. The true 
value probably lies between that defined Eqs. 149 and 102. 

Mr. Krynine’s reference very deep ground-water stream reveals keen 
concept the mathematical premises necessary make solution possible. 
The depth the impervious stratum and its slope may never have been ex- 
plored. fact, clear-cut, impervious stratum may not even exist. In- 
stead, the aquifer may terminate merged zone which, under the influence 
the critical surcharge, has broken down into quasi-impervious soil. 
effective, impervious stratum must then assumed, preferably the form 
horizontal plane. must not overlooked, however, that the pressure con- 
ditions, even deep ground-water stream, open aquifer, may de- 
cisively affected through the creation spreading mound commensurate 
height. shown Figs. and 23, such stream, regardless how mighty, 
may checked strip spreading. The original ground-water flow down- 
stream from the spreading ground then replaced infiltrated water from the 
latter during the first, and part the second, phase. axiomatic that 
accomplish this replacement the infiltrated flow per unit length spreading 
ground must greater than the corresponding initial ground-water flow, which 
reveals the fact that the spreading ground itself control even the mound 
never rises the surface. Hence, spreading grounds are powerful hydraulic 
and must operated with proper appreciation that fact. 

The progressive transmission motion the spreading wave from its 
surface the initial ground-water stream any depth through the medium 
shearing stresses telling example Mr. Krynine’s insight into these phe- 
nomena. His reference the fruitfulness the parabolic bottom used the 
model doubt closely related the significance shearing stresses. 
this connection well realize, however, that the mathematical treatment, 
based the Boussinesq equation general and Eq. (heat flow analogue) 
particular afflicted with physical inconsistency. This due the 
neglect inertia forces. implies instantaneous instead progressive ac- 
celeration through the spreading mound initial ground-water bodies and 
instantaneous transmission pressure the control infinity, during the 
first phase. The introduction the intuitive, unique flow function Part 
Section tends mitigate the aforementioned inconsistency that finite 
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wave length associated with finite time element. interesting note 
from Fig. that equations based the latter concept led closer agreement 
with the test results than those based Eq. 


conclusion, the writer wishes express his sincere appreciation the 


discussers the time and effort which were necessary compose their valuable 
contributions this paper. the light the combined efforts hoped 
that ground-water phenomena related spreading will more readily under- 
stood the future than they have been the past. 
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